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Abstract 

A negative differential resistance (NDR) in nanotransport is often ascribed to electron correla- 
tions. We present a simple example revealing that finite electrode bandwidths and energy depen- 
dent electrode density of states can cause a significant NDR, which may occur even in uncorrelated 
systems. So, special care is needed in assessing the role of electron correlations in the NDR. 



PACS numbers: 73.63.-b, 85.35.Be, 85.35. Gv, 85.65. +h 



The fact that the current-voltage (/-V) characteristics of the dc-transport can exhibit a 
negative differential resistance (NDR) in systems described within a single-particle picture, 
and is not necessarily related to electron correlations is well known in semiconductor physics. 1 
However, in the nanophysics community the NDR in the I-V curve is often ascribed to 
(presumably strong) electron correlations. In fact, some calculations performed on simple 
but nontrivial models of correlated electrons, like the interacting resonant level model, found 
no NDR effect far away from resonance,-^ while other calculations revealed a more^ or less^ 
pronounced NDR effect at resonance. At the end of this note, we shall return to the NDR 
effect within the interacting resonant model. Beforehand — and this is the main aim of 
the present work — we want to emphasize that other, more common sources of the NDR 
are relevant for nanotransport as well. Therefore, special care is needed if one attempts to 
ascribe the NDR to electron correlations. 

The naive "argument" behind the confusion that the NDR is an electron correlation effect 
seems to be the following. Within the Landauer approach of the transport in uncorrelated 
systems, the current resulting from the imbalance between the source and drain chemical 
potentials us = £f + and /io = Ef — eV s d/2 is expressed as an integral of the 

transmission coefficient T(e) over energies from e = fin to e = fi$- An NDR cannot occur 
because the current monotonically increases, since the integrand is positive [T(e) > 0] and 
the integration range increases as the voltage V s d becomes higher. 

To illustrate that this is not the case, let us consider a two-terminal setup (Fig. [1]), 
consisting of a nanosystem [quantum dot(s) or molecule(s)] linked to semi-infinite leads 
(source and drain) at zero temperature. For simplicity, their bandwidth At as well as their 
coupling to (say) the dot r will be supposed to be identical. By gradually rising the source- 
drain voltage V s d starting from V s d = 0, the drain current I s d will first progressively increase 
because the energy window AE of the (elastic) electron tunneling processes allowed by 
Pauli's principle becomes broader (Fig. [T^,). However, further increasing V s d beyond half of 
the electrode bandwidth (eV* d = 2t) will diminish this energy window (Fig. [lb), and this will 
be accompanied by a current reduction, which becomes more and more pronounced as the 
electrode band edge is approach. For eV s d > At, elastic tunneling is no longer possible, and 
the current is completely blocked (I s d = 0). This fact that the current I s d should diminish as 
V s d exceeds V* d and is completely suppressed above the band edge (At) applies for a general 
two-terminal setup for a sufficiently weak hybridization r = 2r 2 /t. 
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FIG. 1: (Color online) Schematical representation of a typical two-terminal setup. By gradually 
increasing the source-drain voltage V s d the energy window AE of the allowed elastic tunneling 
processes (a) increases for eV S( i < 2t, but (b) beyond the point eV s d = 2t (electrode half-bandwidth) 
it decreases. Elastic tunneling cannot occur for eV s d > At. 

To make the analysis more specific, let us consider a point contact (noninteracting reso- 
nant level) model, wherein the nanosystem consists of a single nondegenerate energy level e g 
linked to one-dimensional semi-infinite electrodes. The second-quantized Hamiltonian reads 

H = -t ^2 (4 c i-i + h - c -) + c l c i 

K-l K-1 

-t (4 c i+i + h -°) + 4 c i (!) 

i>i i>i 

+e g clc - r (cLiO) + c\c + h.c)j . 

As usual, we set t = 1 and Ep = 0. We assume e g > (n-type conduction) for simplicity, 
but because model possesses particle-hole symmetry, one can replace e g by \e g \ below. 
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The electrode-dot coupling r yields well known expressions of the embedding self-energies 
E x (e) = A x (e) - iT x (e)/2 (x = S,D), wherein 



A x (e) = A(e - fi x ); T x (e) = T(e - /i x ) 



r 2 e „, , r 2 



A(e) = w ;T(e) = ¥ VW^e(2t-\e\). (2) 

They can be inserted into the Dyson equation 

G - 1 (e) = e-e 9 -i:s(e)-^ D (e) (3) 

to obtain the retarded Green function G(e) of the embedded dot. With the aid of the latter, 
the electric current can be expressed as (electron spin is disregarded) 



l sd — r 
n 



f m deT{e) = \ r deT s (e)T D (e)\G(e)\ 

Jrin tl J/J-n 
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e T D {e)T s {e 
as 



h 



l-'D 



e-e ff -A(e))l 2 + r(e)V4 



(4) 



where Y{e) = T D (e) + T s (e) and A(e) = A D (e) + A s {e). 

I-V characteristics computed exactly by means of Eq. (Tj0) at resonance (e g = 0) are de- 
picted by the thick lines in Fig. |2J These curves show that, indeed, the current is suppressed 
as the bias approaches the bandwidth and disappears beyond eV s d > 4t. Away from res- 
onance (s g 7^ 0), a new aspect is visible in Fig. [3^. The current vanishes even below the 
bandwidth At. Practically, the suppression is complete at V sc i = At — e g \ beyond this value, 
the I-V curves only exhibit negligible tails of widths ~ Tq = 2r 2 /t. On the other side, 
the exact I-V characteristics of Figs. [2] and [3|i reveal that the current decreases well before 
reaching the value V* d = 2t/e, which one could expect from Fig. [TJ This demonstrates that 
the finite bandwidth effect discussed above is only one reason why the NDR should occur. 
Significant physical insight can be gained by examining three limits of Eq. (j3J): 
(i) One can approximate the embedding energies by their values at £ = \i x (£s,d — 
— iTq/2) in the whole integration range, which means to simply ignore the 9 step functions 
in Eq. (j2J). One then gets the current 

run-^o (^eV-2e g , eV + 2e g 



I 7 = _Ji arctan 9 - + arctan ■ 9 - . (5) 

sd h V 2r 2T J K ' 

As this amounts to assume that the electrode bandwidth is the largest energy scale (more 
precisely, for V s d, £ g , r<t), Eq. (jSJ) is usually referred to as the wide band limit. 
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FIG. 2: (Color online) I-V curves at resonance {e g = ef = 0) for r = 0.05; 0.1; 0.2; 0.4 computed 
exactly (thick lines) and within approximation (i) described in the text (thin lines). Current I s d 
in units I s sd = neTo/h. 

(ii) Next, one can compute the current using the electrode density of states (DOS) T x for 
e = fi x , but unlike above, considering the Heaviside 9 functions in Eq. (J2J) 



where A± = [min(eV^, At — eV s d) ± 2e g ] x (1 — r 2 /t 2 ). Similar to approximation (i), the 
electrode DOS is assumed constant, but the fact that the electrode bandwidths are finite 
(the main physical aspect underlying Fig. [T|) is taken into account by this approximation. 

(hi) Because the main contribution to the integral in Eq. (j3J) comes from the pole of the 
Green function of the isolated dot, one can use the embedding energies calculated at e = e g . 
In fact, this approximation yields very accurate I-V curves, which are not shown because 
they could be hardly distinguished from the exact curves within the drawing accuracy of 
Figs. El EK, HI and [5j More instructive is however to furthermore assume that the voltage 
V s d is sufficiently high and extend the integration in Eq. (j?J) from — oo to +oo. The result is 



I-V curves in the limit (i) are depicted in Figs. |2] (thin lines), [3]d, and HJ They show a 
monotonically increasing current, which exhibits a step at eV s d — 2e g of width SV s d increasing 
with r and rapidly saturates at an ^-independent value J| d = iieT /h. Such curves are 





high 
sd 



e T(e g -eV/2)T(e g + eV/2) 
hT(e g - eV/2) + T(e g + eV/2)' 



(7) 
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FIG. 3: (Color online) I-V curves out of resonance for r = 0.1 (Tq = 0.02) computed (a) exactly and 
(b) within approximation (i) described in the text for e g = 0; 0.2; 0.4; 0.6; 0.8; 1 (values increasing 
downwards). Current in units et/h. 

usually shown in textbooks, and this feeds the lore of the absent NDR in uncorrelated 
systems. 

What is wrong with the naive argument against the NDR in uncorrelated systems is that 
the transmission is not independent of V s d- The V^-dependence enters via the electrode 
densities of states T SD [cf. Eq. fl2])]. 

On one side, this dependence is considered by the 9 functions of Eq. (j2J), which diminish 
the window of allowed tunneling processes. Approximation (ii) that accounts for this yields 
two qualitatively correct results: an NDR beyond V^, where the predicted I-V curve exhibits 
a cusp (Fig. EJ) and a vanishing current for eV s d > At. Quantitatively, the NDR onset (at 
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Vsd = V* d ) is unsatisfactory; compare these approximate curves (label fb) with the exact 
ones in Figs. H] and El The NDR occurs well below the point predicted by this approximation. 

On the other side, not only the 9 functions, but also the e-dependence of the electrode 
DOS [the square roots in Eq. (|2J)] is important. It is this fact that makes the finite bandwidth 
argument incomplete. The e-dependence of T s,d is accounted for within approximation (iii). 
The comparison with the exact curves (Fig. H]) reveals an excellent agreement at sufficiently 
higher voltages (as assumed within this approximation) and demonstrates that, to describe 
quantitatively the NDR, one has to consider both the allowed energy window, which is finite, 
and the energy dependence of the electrode DOS. 

In Fig. HI we present exact I-V characteristics from Eq. (J4J) along with those computed 
within the three aforementioned approximations, Eqs. (jSJ), flS]), arid (j7|). As visible there, 
approximation (i) is accurate for lower voltages, while approximation (iii) is accurate for 
higher voltages. The crossover occurs at a voltage , which can be identified with the 

NDR onset. This value can be obtained by equating 

i i 7(v s r R )=iT(v s N d DR )- (8) 

Curves for V^ DR are presented in Fig. |6j They show that for situations not very far away 
from resonance and sufficiently weak electrode-dot couplings r, V S ^ DR is considerably smaller 
than the value eV* d = 2t expected from the finite bandwidth argument. The significant 
departure of the NDR onset predicted exactly and within approximation (ii) is also clearly 
depicted in Fig. [51 For smaller r's one can deduce an analytical estimate (c ~ 4) 

V s N d DR ^2e g + c(tT*)V 3 . (9) 



Interesting for nanotransport are the electron level(s) not too misaligned with electrode's 
Fermi level; otherwise, as illustrated by the curve for e g = t in Fig. [3^, the current is very 
small. Therefore, the results on V S ^ DR expressed by Eq. (Q and Fig. [6] are perhaps the most 
relevant ones from an experimental perspective. At resonance and realistic parameters (t ~ 
1 eV, r 1 meV-), Eq. yields V^ DR ~ 40 meV. Based on this estimate, we argue that the 
NDR discussed here can be observed. On one side, correlations are important only at much 
lower voltages; in single-electron transistors,- the relevant scale is the Kondo temperature Tk 
[eV s d ~ k B T K < 0.1 meV). For voltages of tens of mV, correlation effects (e. g., Kondo's) are 
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FIG. 4: (Color online) 7-V curves for r = 0.1 computed exactly and within the approximations 
described in the text: (a) at resonance Ed = and (b) out of resonance, = 0.2. Current in units 
et/h. Labels as in Eqs. (|5|), ([6]), and (J7]). 



supprressed; the present uncorrelated limit is justifiable. On the other side, the estimated 
NDR onset voltages (~ 10 mV) are much lower than the electrode bandwidth (~ leV), 
and a material damage prior to the NDR onset can be ruled out. For Si-based SETs, the 
material can support even much higher values, V s d ~ IV. 10 So, we hope that the present 
estimate will stimulate experimentalists to search NDR effects at moderate V s d- Again quite 
relevant for experiments, the NDR onset can be controlled by tuning the level's energy e g 
with the aid of a gate potential. Gating methods were routinely employed for nanosystems 
in the past£ and recently also in molecular transport.- 1 In (weakly-correlated) molecules, 
the level e g would be either the highest occupied molecular orbital (HOMO)ii or the lowest 
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FIG. 5: (Color online) I-V curves on resonance (e g = 0) for the three electrode-dot couplings 
r specified in the inset computed exactly and within approximation (ii) described in the text 
(label fb). Notice that the latter exhibit a cusp at eV s d = 2t that marks the NDR onset in this 
approximation, which can be substantially higher than the exact NDR onset. 




FIG. 6: (Color online) Curves for the NDR onset voltage V^ DR computed from Eq. ([8]) for several 
level energies e g . Notice that for smaller electrode-dot couplings r and not too far away from 
resonance, V^ DR is significantly smaller than 2 (half of electrode's bandwidth). 

unoccupied molecular orbital (LUMO, as in Fig. [1]), depending on which is closer to Sp. 
There, r ~ leV and \e g \ ~ leV.— So, the NDR-onset [cf. Eq. (JSJ) and Fig. [6] is expected at 
V s d- values of a few eV, slightly higher than used in experiment.— 

The present analysis can be extended without difficulty to nanosystems/molecules with 
several "active" electron levels. As long as these levels £ g2 , ■ ■ ■ are well separated energet- 
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ically and the hybridization is weak enough (a different situation can also be encountered, 



see Ref. 



12l ). they manifest themselves as current steps at the voltages eV s d ~ 2e g i, 2e g 2, 



However, even in this case the finite electrode bandwidth and the energy dependence of the 
electrode DOS remain possible important sources of an NDR. 

Similar to other situations encountered in nanotransport,— ^ we believe that the results 
for uncorrelated systems are instructive and could also be useful to correctly interpret the 
nanotransport in correlated systems. In the present concrete case, they could help to unravel 
the physical origin of the NDR. In the light of the present analysis, it is plausible to ascribe an 
NDR as an electron correlation effect in cases where the NDR was found within calculations 
to a correlated nanosystem carried out within the wide band limit. This is, e. g., the case of 
Refs-Oandlsl, where a weaker NDR effect was obtained at resonance at stronger Coulomb 
contact interactions. As suggested by Fig. [3l the farther away from resonance, the more 
is the NDR onset pushed towards higher voltages (eV s ^ DR > 2|e fl |). The values of V s d 
chosen in the figures shown in Ref. [3 do not belong to this range and the absence of an 
NDR could be related to this fact. Unlike the wide (infinite) band limit assumed in the 
aforementioned references, a discrete model of the electrodes, with a finite bandwidth At, 
exactly as in Eq. (JT]), has been utilized for the numerical calculations of Ref. |4] at resonance. 
The I-V curves reported there exhibit a pronounced NDR effect. However, in view of the 
finite bandwidth assumed in that work, attributing this effect to electron correlations at 
rather high voltages should be made with special care. We believe that in order to interpret 
this effect reliably, one should first carefully subtract the contribution to the NDR due to 
the finite bandwidth and the energy dependent electrode DOS discussed above. 
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fully acknowledged. 
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